The spherical collapse model in time varying vacuum cosmologies 
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We investigate the virialization of cosmic structures in the framework of fiat FLRW cosmological 
models, in which the vacuum energy density evolves with time, fn particular, our analysis focuses 
on the study of spherical matter perturbations, as they decouple from the background expansion, 
"turn around" and finally collapse. We generalize the spherical collapse model in the case when the 
vacuum energy is a running function of the Hubble rate, A = A(H). A particularly well motivated 
model of this type is the so-called quantum field vacuum, in which A(H) is a quadratic function, 
A(H) — no + ni H 2 , with no 7^ 0. This model was previously studied by our team using the latest 
high quality cosmological data to constrain its free parameters, as well as the predicted cluster 
formation rate. It turns out that the corresponding Hubble expansion history resembles that of 
the traditional ACDM cosmology. We use this A(i)CDM framework to illustrate the fact that the 
properties of the spherical collapse model (virial density, collapse factor, etc.) depend on the choice 
of the considered vacuum energy (homogeneous or clustered). In particular, if the distribution of 
the vacuum energy is clustered, then, under specific conditions, we can produce more concentrated 
structures with respect to the homogeneous vacuum energy case. 
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INTRODUCTION 



Several cosmological observations (supernovae type la, 
CMB, galaxy clustering, etc.) have converged to a 
paradigm of a cosmic expansion history that involves a 
spatially flat geometry and a recently initiated acceler- 
ated expansion of the universe (see [l|-|6[ and references 
therein). From a theoretical point of view, an easy way to 
explain such expansion is to consider an additional energy 
component, usually called dark energy (DE) with nega- 
tive pressure, that dominates the universe at late times. 
The simplest DE candidate corresponds to the cosmolog- 
ical constant (see [7H9| f° r reviews). An elegant model 
that fits accurately the current observational data, is the 
spatially flat concordance ACDM model, which includes 
cold dark matter (DM) and a cosmological constant, A. 
However, the A model suffers, among other [lOJ . form 
two fundamental problems: (a) The "old" cosmological 
constant problem (or fine tuning problem) i.e., the fact 
that the observed value of the vacuum energy density 
(pA = Ac 2 /8ttG ~ 10 -47 GeV 4 ) is many orders of mag- 
nitude below the value found using quantum field theory 



ll| i.e., the fact that 



[7|, and (b) the coincidence problem 
the matter energy density and the vacuum energy density 
are of the same order (just prior to the present epoch), 
despite the fact that the former is a rapidly decreasing 
function of time while the latter is stationary. 

Many authors have attempted to solve the above prob- 
lems (see |a |9|, Il2[ and references therein) , the key ap- 
proach being the replacement of the constant vacuum 
energy either with a time evolving DE (quintessence and 



the like [8] ) , or alternatively with a time varying vacuum 
energy density, p\(t) [13H18J . In the original scalar field 
models [l9J] and later in the quintessence context, one can 
ad-hoc introduce an adjusting or tracker scalar field <f> [20] 
(different from the usual SM Higgs field), rolling down 
the potential energy V(<p), which could mimic the DE 
[a IqL l2l| - [24J . However, it was realized that the idea of a 
scalar field rolling down some suitable potential does not 
really solve the problem because <f> has to be some high 
energy field of a Grand Unified Theory (GUT), and this 
leads to an unnaturally small value of its mass, namely 
one which is beyond all conceivable standards in Parti- 
cle Physics. As an example, utilizing the simplest form 
for the potential of the scalar field, V{4>) — m 2 4> 2 /2, the 
present value of the associated vacuum energy density is 
Pa = (V((f>)) ~ 10~ n eV 4 , so for (0) of order of a typical 
GUT scale near the Planck mass, Mp ~ 10 19 GeV, the 
corresponding mass of <\> is expected in the ballpark of 
m ~ # ~ 1CT 33 eV. 

Notice that the presence of such a tiny mass scale in 
scalar field models of DE is generally expected also on the 
basis of structure formation arguments |25l |26| ; namely 
from the fact that the DE perturbations seem to play 
an insignificant role in structure fornration for scales well 
below the sound horizon. The main reason for this homo- 
geneity of the DE is the flatness of the potential, which 
is necessary to produce a cosmic acceleration. Being the 
mass associated to the scalar field fluctuation, propor- 
tional to the second derivative of the potential itself, it 
follows that ?Ti0 will be very small and one expects that 
the magnitude of DE fluctuations induced by <f> should be 
appreciable only on length scales of the order of the hori- 



zon. Thus, equating the spatial scale of these fluctuations 
to the Compton wavelength of 4> (hence to the inverse of 
its mass) it follows once more that m^ < Hq ~ 10 -33 eV. 
All in all, it appears that the problem that one is creating 
along with the introduction of <f> is far more worrisome 
than the problem one is intending to solve, for one is 
postulating a mass scale which is 30 orders of magnitude 
below the mass scale associated to the value of the vac- 
uum energy density (ttt-a = p A ' ~ 2.3 x 10 -3 eV). 

The analysis of the recent cosmological observations 
indicates that the DE equation of state (EOS) parame- 
ter w(= Pde/pde) is close to —1 to within ±10%, if it 
is assumed to be constant [Iref, whilst it is much more 
poorly constrained if it varies with time [3] . More than 
two decades ago, Ozer & Taha [27( proposed a time vary- 
ing A as a possible candidate to solve the two fundamen- 
tal cosmological puzzles, see also [28|, [29[ and references 
therein. In this cosmological paradigm, the dark energy 
EOS parameter w is strictly equal to -1, but the vac- 
uum energy density (or A) does evolve with time. Of 
course, the weak point in this approach is the unknown 
functional form of A(i), which is however also the case for 
the vast majority of the DE models. Indeed, in the afore- 
mentioned A(i) models, the evolution law is purely phe- 
nomenological 30] , without a concrete link to fundamen- 
tal physics, such as the Quantum Field Theory (QFT) in 
Curved Space time |3l| . As emphasized in [32| , a com- 
pletely consistent formulation along these lines should 
eventually be developed, and such investigations could 
well be at the heart of one the most important endeavors 
of theoretical cosmology in the years to come. Therefore, 
the study of cosmic perturbations in these models is very 
important 33, 34] as they might reveal surprises not fore- 
seen in the context of the scalar field models. The new 
effects may have impact both on the cosmological and the 
astrophysical domains. While we recently analyzed the 
potential implications on the former [35j , here we focus 
on the latter domain. 

A pioneering QFT fundamental approach to variable 
A models was actually proposed long ago within the 
context of the renormalization group (hereafter RG) in 
|36l 1371 ] . Later on, the RG-running framework was fur- 
ther explored from different points of view in [L3l - [l5| . 
and a more systematic presentation from the viewpoint 
of QFT in curved space-time by employing the standard 
perturbation RG-techniques of Particle Physics appeared 
in [l4j, |l6j. Subsequent elaborations, and comparison 
with the observational tests, confirmed the phenomeno- 
logical viability of this approach [38|-[44| . 

In the class of RG models we shall focus on, the vacuum 
energy density is expected to vary with time according to 
the law [Hl3llJ4l"j]: A = n +n 2 H 2 (hereafter called the 
Aug model or quantum field vacuum model). As already 
mentioned, in Ref. [35J we have investigated thoroughly 
the global dynamics of this cosmological model (together 
with various alternative A(i) models), in the light of the 
most recent cosmological data. However, a serious issue 
here is how the main features of the largest collapsed 



cosmic structures, i.e., galaxy clusters, are affected by a 
running vacuum energy density. We have argued above 
that this problem can be addressed in scalar field mod- 
els of the DE, but only at the expense of admitting ex- 
tremely tiny mass scales which are uncommon in Particle 
Physics. In this paper, we wish to further explore the al- 
ternative option in which the DE component is a time 
evolving cosmological term A = A(i), and in this way to 
assess if the clustering properties of the vacuum energy 
can shed some light on the fundamental issue of structure 
formation. 

The so called spherical collapse model [45] , which has 
a long history in cosmology, is a simple but still a funda- 
mental tool used to describe the growth of bound systems 
in the universe via gravitation instability [461 ] . In the last 
decade many authors have studied the small scale dynam- 
ics using this model and found that the main features of 
the cosmic structures (collapse factor, virial density, etc) 
can potentially be affected by the presence of dark energy 
[25J 12a . [47H66l | . The aim of the present work is to general- 
ize the spherical collapse model within the variable Aug 
cosmological model, in order to understand non-linear 
structure formation in such cosmologies and investigate 
the differences with the respect to the expectations of the 
concordance ACDM cosmology. 

The structure of the paper is as follows. The basic the- 
oretical elements of the problem are presented in section 
2, where we introduce [for a spatially flat Friedmann- 
Lemaitre-Robertson- Walker (FLRW) geometry] the ba- 
sic cosmological equations. In section 3 we generalize the 
virial theorem in the case of the QFT A(t) cosmologi- 
cal model. Section 4 outlines the theoretical analysis of 
the spherical collapse model in which the vacuum energy 
density varies with the cosmic time, and in section 5 we 
compare the corresponding theoretical predictions of the 
different models and present a first attempt to use ob- 
servational data to constrain the different models. We 
draw our conclusions in section 6. In the appendix A 
we remind the reader of some basic elements of the con- 
cordance ACDM model in order to appreciate the fact 
that the Arc cosmology is an interesting extension of 
the standard model. Finally, in appendix B we provide 
some basic mathematical formulae, while in appendix C 
we provide accurate fitting formulae for a few important 
parameters, ie., the density contrast at the turn around 
redshift and at the epoch of virialization, which do not 
have a simple fully analytical form. 



2. COSMOLOGY WITH A TIME DEPENDENT 
VACUUM 

The cosmological constant contribution to the curva- 
ture of space-time is represented by the A g^ term on the 
l.h.s. of Einstein's equations. The latter can be absorbed 
on the r.h.s. of these equations 
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Rfiv — t,9^vR — 87rG T^v , 



(1) 



where the modified T^ v is given by T M „ = T^ v + g^ v p A - 
Here pa = A/(8ttG) is the vacuum energy density associ- 
ated to the presence of A (with pressure Pa = ~Pa), 
and T^ v is the ordinary energy-momentum tensor of 
isotropic matter and radiation. Modeling the expand- 
ing universe as a perfect fluid with velocity 4- vector field 
U M , we have T^ = -P m g M „ + (p m + P m )U M U,„ where 
p m is the proper isotropic pressure of matter-radiation 
and P m is the corresponding pressure. Clearly the mod- 
ified T^ v defined above takes the same form as T^ v with 
p_ = Pm + PA and p = P m + P A = P m - p A , that is 
T^u = -pg^v + (p + p)U^U u . Explicitly, 



Tfiv = (p\ - Pm) 9^u + (Pm + Pm)UfJJ v 



(2) 



With this generalized energy-momentum tensor, and in 
the spatially flat FLRW metric ds 2 = dt 2 — a 2 (t)dx 2 , the 
gravitational field equations boil down to Friedmann's 
equation 



H 2 



InG 8ttG . . 

-5— P = ~^- {Pm + PA) , (3) 



and the dynamical field equation for the scale factor: 

a AttG , , A-kG , „ , . 

- = —(p + 3p) = — (p m +3P m -2p A . (4) 

a 3 3 

Let us next contemplate the possibility that pa = Aa(£) 
is a function of the cosmic time. This is allowed by the 
Cosmological Principle embodied in the FLRW metric. 
The Bianchi identities (which ensure the covariance of 
the theory) then imply v M T^v = 0. With the help of 
the FLRW metric, the previous identity amounts to the 
following generalized local conservation law: 



Pn 



p\ + 3H{p m + P m + pa + Pa) = . 



(5) 



where the over-dot denotes derivative with respect to the 
cosmic time. The above equation can also be derived by 
combining eqs. ([3]) and (0| since it is a first integral of 
the equations of motion. Notice that we keep G strictly 
constant, and therefore the assumption p\ 7^ necessar- 
ily requires some energy exchange between matter and 
vacuum, e.g. through vacuum decay into matter, or vice 
versa 1 . This possibility was first considered by M. Bron- 
stein in a rather early paper [68] . 

Let us remark that the EOS of the vacuum energy 
density maintains the usual form Pa(£) = —p\(t) — 
— A(£)/87rG despite the fact that A(t) evolves with time. 
In the matter dominated epoch (P m — 0) , eq. (JS|) leads to 



1 There exists also the possibility that the vacuum is time evol- 
ving and nevertheless non-interacting with matter. In this case, 
however, either the DE has another component apart from A— 
see the AXCDM framework of Il7l — or Newton's coupling is also 
time-varying, i.e. G ± H, IMlM E3] . 



the following energy exchanging balance between matter 
and vacuum: 



3Hp n 



-PA- 



(6) 



The second Friedmann's equation (j4]) is formally un- 
changed by the presence of a time-variable vacuum en- 
ergy, and in the matter epoch simply reads 



4ttG 



{pm - 2pa) • 



(7) 



At this point it is worth noticing that the effect of having 
a variable cosmological term pa — pA{t) cannot, in gen- 
eral, be described by the simple parameterizations usu- 
ally employed for the effective EOS w — w(t) of the DE, 
in which w depends on two parameters (wq,wi), that 
can constrained using observations [2l|, |69| . The effec- 
tive EOS of a variable vacuum model is in general more 
complicated. This is shown in detail, with specific exam- 
ples, in reference [401 ] . In particular, the vacuum models 
that we are going to consider cannot be described with 
these simple parameterizations. Therefore, the variable 
vacuum models must be studied on their own and con- 
stitute an independent class of DE models. 
Combining equations ([3]) and ([6]), we find: 



H + ±H 2 



A 
4ttGpa = - . 



(8) 



If the vacuum term is negligible, A — > 0, then the solu- 
tion of eq. (jHJ) reduces to that of the Einstein-de Sitter 
model, H(t) = 2/3t, as it should. Similarly, the tra- 
ditional A =const. cosmology (or ACDM concordance 
model) also follows directly by integrating eq. (JU) [see ap- 
pendix A]. Finally, this same equation is also valid for 
A = A(£) when matter and vacuum become coupled, and 
in this case a supplementary equation for the time evo- 
lution of A is needed in order to unveil the dynamics of 
this model. It is interesting to mention here that the 
link in eq.© between p m and A is important because 
interactions between DM and DE could provide possible 
solutions to the cosmological coincidence problem. This 
is the reason for which several papers have been pub- 
lished recently in this area [13, Lis IZO] ; proposing either 
that the DE has various interacting components or that 
the DE and DM could be coupled. In the following sub- 
sections, we briefly introduce (for more details see [351 ]) 
the cosmological models used in this study. 



2.1. The A(i) model from quantum field theory 

In this scenario, we use the vacuum solution proposed 
in [IT], |38h41| ] using the renormalization group (RG) in 
quantum field theory (hereafter A#g model). The model 
is characterized by the evolution law: A(H) = no+ri2 H 2 , 
in which both coefficients no and 712 are non- vanishing. 
An equivalent and more convenient parametrization is 



A(H) = A + 3 7 (H 2 



Hi 



), 



(9) 



where 7 is a constant, which can be positive or negative 
but small: I7I < l/Ylir 35]. It determines the amount 
of running of A(t). Note that the form (|9|) is convenient 
because the vacuum energy density is normalized to the 
present value: Ao = A(Hq) — 3£IaHq. As we will see, all 
the equations and conditions derived below are equiva- 
lent to those of the concordance ACDM cosmology (cf. 
appendix A) and reduce exactly to them for 7 = 0. 

From equations flSJ and (HJ) we can easily derive the 
corresponding Hubble flow as a function of time [35[ : 



H(t)=H Q J— — 1 coth 
V !-7 



-iW(^A-7)(l-7)* 



The scale factor of the universe, a(i), evolves as 
"3 



a(t) = cii sinh 3 < 1 -T) 
where 

«i 



■*W(fiA-7)(l-7)* 



/ n ri 



\Qa -7 



(10) 



, (11) 



(12) 



Inverting eq. (|11[) we determine the cosmic time as a func- 
tion of the scale factor: 



t(a) 



3^ /2 (l-7)if 



sinh 




A 3(l- 7 )/2 



where we have introduced 

1 L r n 



- £ r, 



1 



'7 



^A 



^A -7 

1-7 



(13) 



(14) 



Let us define the Hubble rate normalized to its current 
value, E(a) = H(a)/H . Eliminating the cosmic time 
from equations (ITU1) and (jTTJ), one can prove that 



Ei2/ \ ^A 7 , "m 

E ( a ) = -j + "J a 

1 — 7 1 — 7 

a -3(l-7) . 

= i + n m 



-3(1-7) 



1 



7 



(15) 



Using p m = aHdp m /da, we can trivially integrate the 
previous equation in the scale factor variable, yielding 



Pm{a) = Pmoa 3(1 7) 



(19) 



where p m o denotes the matter density at the present 
time (a = 1), and therefore O m = Pmo/Pco, where 
p c x) = 3H 2 /(8irG) is the current critical density. The 
previous equation can also be rewritten by considering 
the instantaneous critical density p c (a) when the scale 
factor is a — a(t), i.e. p c (a) = 3H 2 (a) / (8irG) . In fact, 
defining fi m (a) = p m {o) / p c {a) it is easy to see, with the 
help of ([T9")) and the definition of E(a), that 



O m (a) = 



il, 



-3(1-7) 



E 2 {a) 



(20) 



Finally, upon inserting (fT9")) in © and integrating once 
more in the scale factor variable, we arrive at the ex- 
plicit expression for the evolution of the vacuum energy 
density: 



A(o) = A + 8ttG 



7 PmO 



1 



7 



a-^ 1 -^ - 1 



(21) 



It is important to emphasize from eq. (|19l) that the matter 
density does no longer evolve as p m (a) = p m Qa 3 , as it 
presents a correction in the exponent. This is due to the 
fact that matter is exchanging energy with the vacuum 
and this is reflected in the corresponding behavior of A(a) 
in eq. (|2T|) . As expected, for 7 — > (fi m ~ fi m , 17a ~ ^a) 
all the above equations reduce to the canonical form 
within the concordance model (for more details, see the 
appendix A). Clearly, the usual A-cosmology is a par- 
ticular solution of the Arc model with 7 strictly equal 
to 0. Throughout the rest of the paper we shall em- 
ploy the statistical results for this model obtained in [35j | 
from a simultaneous fit to the latest SNIa+BAO+CMB 
data, namely Q m = 0.28;to - oi (or O m ~ 0.281) and 
7 = 0.002 ±0.001. 



GENERALIZATION OF THE VIRIAL 
THEOREM 



E 2 (a) = tt t 



•a. 



-3(1-7) 



Notice that both sets of cosmological parameters sat- 
isfy the standard cosmic sum rule: 



- ';, 



Q\ = 1 = Vi r 



£1/ 



(17) 



Let us also consider the evolution of the matter and 
vacuum energy densities in this model. Starting from 
the conservation law [see eq(5] and utilizing eqs.([3]) and 
(0), we arrive at a simple differential equation for the 
matter density, 



p, n + 3Hp m = 3jHp„ 



Recall that for systems with potential energy of the 
(16) form U oc R n the contribution to the virial condition 
is 2T — nil = 0, where T is the kinetic energy. For 
gravity, n = — 1, whereas for constant vacuum energy 
n = 2. Indeed, from Newton's limit of Einstein's equa- 
tions in the presence of a A term, i.e. Poisson equation 
V 2 $ = AirG(p m — 2p\), it follows that the potential as- 
sociated to a constant A is ^a = —(1/6) AR 2 . Thus, the 
nominal virial condition when there is a constant vacuum 
energy reads 2T + Uq — 2U\ = 0, where U~g and Ua are 
the potential energy and the vacuum potential energy, 
respectively, for an isolated system. However, for non- 
constant A this recipe may not hold anymore. Therefore, 
(18) the link between the time-dependent vacuum and matter 



is expected to modify non-trivially the previous form of 
the virial theorem. 

In this section, we generalize the virial theorem by ta- 
king into account the presence of the coupling between 
the vacuum and the matter energy densities, which leads 
to a variable A(t). In particular, we have to modify the 
well know Layzer-Irvine equation, which describes the 
flow to virialization [46[. As we have already stated in 
section 2.1 (see eqj6]), the matter is exchanging energy 
with the vacuum and this is reflected in the correspond- 
ing matter continuity equation (|18[) . Of course, the con- 
tinuity equation holds for the total cosmic fluid. Since, 
we are interested for cosmic structures which live in high 
density environments, it is fair to consider that the cor- 
responding inhomogeneous density n m is far from the 
background homogeneous density, p m . 

The total velocity of the fluid elements, v, is given by: 
V • (n m v) = 3Hn m + V • (n m u), where u is the peculiar 
velocity. Note that v is the sum of the global Hubble 
expansion plus the line of sight peculiar velocity, u. The 
continuity equation for n m , in the presence of a running 
vacuum energy (7 ^ 0), can be now modified as follows: 



+ 3Hn m + V ■ (n m u) = i^Hn^ 



(22) 



Notice that if the background would not be expanding 
(H = 0), then eq. (|2"2"|) boils down to the standard 
non-relativistic continuity equation for a fluid of parti- 
cles flowing with a velocity distribution u. Based on a 
Newtonian formulation, the acceleration due to gravity 
in comoving coordinates is 



d(au) 
dt 



-aV$ 



(23) 



where $ is the gravitational potential. Multiplying 
eg. ([23]) by the quantity an m u and using the continuity 
cq. ([22l) together with the overall Poisson equation [46| 
V 2 $ = 4:TrG(n m — 2p\), we can integrate over the volume 
in order to define the generalized Layzer-Irvine equation. 
After some calculations we arrive at 



dt ' v dt 

where U = U G - 2U A , 



2 



HU) + 6ja 2 HU (24) 



and 



Uo = --G 



U A = --G 



u 2 n m dV , 

n m (x)n m (x) 
\x — x'\ 

n m (x)pA(x') 
\x — x'\ 



dVdV 



dVdV 



(25) 



(26) 



(27) 



For a spherical mass fluctuation M ~ 47rn m i? 3 /3, one 
can show that the above potential energies become 



16tt 2 G f R 4 2 . . , 

U G = — / x n m (x)dx 



3GM 2 
5R 



(28) 



and 



U A = 



16tt 2 G f R 
3~~ 



AM n 
x pA{x)n m {x)dx = — — -R (29) 



where the last equality holds for a homogeneous vacuum 
energy A = A (a) [see [56|, [58| and section 5]. 

Now, for a system that reaches the equilibrium (virial 
regime, and 14 = T = 0) we derive the condition 



(2 - 3 7 )T + (1 - 67) (U G ~ 2U A ) = 



(30) 



Ug - 2Ua 

T 



1 - 3 7 /2 
I-67 



-2-9 7 + 0( 7 2 ) (31) 



where the last equality is valid for small values of 7. For 
7 strictly equal to zero we obviously recover the nominal 
virial theorem in the concordance cosmology ( c/g ^ 2C/a = 
—2) at it should. In the case of 7 ~ 1/127T [35|], which cor- 
responds to receiving quantum effects on the A running 
from fields just at the Planck scale [lj, |3J|, the above 
ratio in the A(£) cosmology deviates by ~ 12% with re- 
spect to that of the A cosmology. In practice, however, 
7 ~ 1/127T is excluded by the latest fit to the combined 
data, which yields a quite smaller value: 7 ~ 0.002.[35j. 
This value corresponds to quantum effects from GUT 
fields just one order of magnitude below Mp. In this 
case, the deviation from the nominal virial condition is 
only ~ 1%. 
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FIG. 1: The turn around redshift as a function of the virial 
redshift, see eq. H34[) . The solid, long dashed and short dashed 
lines represent the Arc model, for 7 = 0.002, 1/127T and 
— 1/12-7T, respectively. The solid points correspond to the 
concordance ACDM cosmology. Inset Panel: The relative 
fractional difference between the three Ana models and the 
concordance model. 



THE SPHERICAL COLLAPSE MODEL 



Despite its simplicity, the spherical collapse model [45 1 
is still a powerful tool for understanding how a small 
spherical patch of homogeneous overdensity forms a 
bound system via gravitation instability - for a review, 
see e.g. [7l|. Technically speaking, the basic dynamical 
cosmological equation (|Tj) is valid both for the entire uni- 
verse and also for a homogeneous spherical perturbation. 
In the last case, we just replace the scale factor a(t) with 
the radius R(t), and we obtain the so called Raychaud- 
huri equation: 



R 
R 



AnG 



[f>n 



2pAs 



(32) 



where fo = fl\/fl m . 

In the main panel of Figure 1, we present the 
turnaround rcdshift, z t = (1 — at)/a t , as a function of the 
virial redshift z c = (1 — a c )/a c for our Aug model with 
7 = 0.002 (continuous line), 7 = 1/12-7T and 7 = — 1/127T 
(long and short dashed lines, respectively). The concor- 
dance A cosmology is indicated by empty points. The 
relative fractional differences between the A^g models 
and the concordance A model are extremely small, a fact 
which can be appreciated in the inset panel of Fig.l. 

As it is evident there is a tight correlation between the 
two redshifts, which for our Arq model (7 = 0.002) it is 
given by: 



Zt 



1.532z c + 0.751 



(34) 



where p ms and p\ s refer to the corresponding values of 
the matter and vacuum energy densities in the spherical 
patch susceptible of ulterior collapse. 

In order to address the issue of how the time vary- 
ing vacuum energy itself affects the gravitationally bound 
systems (clusters of galaxies), one has to deal in general 
with the following three distinct scenarios, which have 
been considered within different approaches in the lit- 
erature [25J, [26, 47-66]: (i) the situation in which the 
vacuum energy remains homogeneous and only the cor- 
responding matter virializes; (ii) the case with clustered 
vacuum energy, but now assuming that only the matter 
virializes; and, finally, (iii) the case with clustered vac- 
uum energy, considering that the whole system virializes 
(both matter and vacuum components). In this paper, 
we are going to focus on scenarios (i) and (iii) within the 
framework of time varying vacuum energy density. 

From now on, we will call at the value of the scale factor 
of the universe where the spherical overdensity reaches its 
maximum expansion (i.e. when R = 0) and a c the scale 
factor when the sphere virializes, implying that a cos- 
mic structure has formed. Similarly, R t and R c stand 
for the corresponding radii of the spherical overdensity, 
the former being the turnaround (or "top hat") value 
at the point of maximum size, and the latter refers to 
the eventual situation when the sphere has already col- 
lapsed and virialized. Note that due to the coupling be- 
tween the time-dependent vacuum and the matter com- 
ponent one would expect that the matter density in the 
spherical region should obey the same power law as the 
background matter p m (a) ex a -3 ^ 1 ^ 7 ^ (see eq[T9|). Thus, 
Pms oc i? -3 ' 1 / ' denotes the matter density in the spheri- 
cal patch. Analogously, the vacuum energy density in the 
spherical region, p\ s will take the form (|21|) with appro- 
priate replacement of the scale factor with R, see further 
below. 

From the theoretical point of view, the time needed for 
a spherical shell to re-collapse is twice the turn-around 
time, tf ~ 2tt, which implies that (see eqfT3|): 



sinh 



r a 



3(1-7) 



2sinh 



3(1- 



■7) 



(33) 



As an example, assuming that galaxy clusters have 
virialized at the present time, z c ~ 0, the turn around 
epoch takes place at z t ~ 0.75 (or a t ~ 0.57). On 
the other hand, considering that clusters have formed 
prior to the epoch of z c ~ 1.6 (a c ~ 0.38), in which 
the most distant cluster has been found 72], the turn 
around epoch is not really affected by the vacuum en- 
ergy component, i.e. z t ~ 3.2 (or at ~ 0.24). This is to 
be expected, due to the fact that at large redshifts mat- 
ter dominates the Hubble expansion. It is worth noting 
that, at high redshifts, the ratio between the scale factors 
approaches the Einstcin-dc Sitter (il m = f2 m = 1) value: 
a c /a t = {l + z t )/(l + z c )=2 2 / 3 . 

Performing the convenient transformations into dimen- 
sionless variables 



a , R 

— and y = — 
a t Rt 



(35) 



the evolution of the scale factor of the background and 
of the perturbation (see eqs(21 [TBI and |3"2"|) are governed 
respectively by the following two equations: 



H t Q r , 



-3(1-7) 
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y3(l- 7 ) 



(36) 



(37) 



where Hfil mt — ^r^-pm.t, Q m ,t is the matter density pa- 
rameter at the turn around epoch (see eq |2"0|) 2 . Note that 
we have parametrized the matter density in the spheri- 
cal region at the turn around time p ms ,t = CPm.t, with 
respect to the background matter density at the same 
epoch p m> f The parameter £ is referred to as the den- 
sity contrast at the turnaround point. In replacing the 



2 We set U m ,t = f2m(at), with O m (o) given by eg. 1201 . 



quantity p„ 
relation: 



in eq. (|32p we have utilized the following 
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-3(1-7) 
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3(1-7) 
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Therefore, using eas. (|14l21l) we find that 
1 
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3(1-7) 



(38) 



(39) 



p^t being the vacuum energy density at the turn around 
epoch and Da,* = 1 — ^m,t is the effective vacuum density 
parameter at the same time (for definition see eq |20l) . 
Inserting now eq. (|39p into ea. (|36|) . we finally obtain 



x 2 



= H 2 n„ ltt [x 



-1+37 



l I{x)] 



where 



PA,t __ ^A 3(1-7) . 7 



Pm.t 



- -i) 



1 



'7 
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3(1-7) 



(40) 



(41) 



(to derive the latter equality we have used eqs.17,19,21). 
Of course for the Einstein-de Sitter case (Cl m — Q m — 
1 and 7 = 0) the solution of the system formed by eq. (j3"o) 
and eq. (|3"T)) reduces to the well known value of the den- 
sity contrast at the turnaround point: ( = (*¥■) , as 
it should. Within this framework, utilizing both the 
virial theorem (see eq l30p and the energy conservation 
(T c + Ug.c + U\, c = Uc.t + Ux,t) at the collapse time 
and at the turn around epoch respectively we derive the 
following useful relation: 



qiUc,c + Q2Ua,c = U G j + U\ tt 



where 



9i(7) = 



l + 3 7 
2-3 7 



92(7) = 



4-157 

2-37 



(42) 



(43) 



4.1. Homogeneous Vacuum Energy 

In this section, we consider that the vacuum energy 
component of the scale of galaxy clusters can be treated 
as being homogeneous: pAs{a) = Pa{o) = A(a)/8irG 
(hereafter the A-rgh model). Therefore, inserting eq. (j3"9")l 
into eq.p7p. we obtain 



y 



H t Q r 



C. 



2-37 



2ryl(x) 



(44) 



The numerical solution for the £ parameter is pro- 
vided by integrating the main system of differential equa- 
tions, (eqsHOl and 03]), using the boundary conditions: 
(dy/dx) — and y = 1 at x = 1. Following the metho- 
dology of [48| and [65| we provide in the appendix C a 
reasonably accurate fitting formula for £, as a function 
of the main cosmological parameters. 
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FIG. 2: The density contrast at the turn around epoch, ( 
as a function of the turn around redshift. The lines repre- 
sent the following cosmological models: (a) Argh (solid), (b) 
Argci (dot line, 7 S = 0.002) and (c) Argc2 (dashed line, 
7 S = —0.002). The points represent the concordance ACDM 
cosmology. 



Using now, the combined equation (|42p for the poten- 
tial energies (see eas l28l29")) 3 , we obtain a cubic equation 
that relates the ratio between the virial (R c ) to the turn- 
around outer radius (Rt), the so called collapse factor 
(X = R C /R t ): 



? 2 (7HA 3 - (2 + n t )X + 2ft (7) = , 



where 
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(45) 



(46) 



(47) 



(48) 



Finally, solving the cubic eq. (j45j) . we calculate the col- 
lapse factor (see appendix B). In the case of 7 = 
(n c = n t = no), the above expressions take the usual 



3 In view of the fact that Ua = —A(a)MR 2 /10, the time depen- 
dence of the vacuum energy density seems to create a problem 
since the total energy of the bound system is not conserved. 
However, one can show that if the value of I7I is less that 0.01 
then the problem of energy conservation does not really affect the 
virialization process and thus eq. ||42|I remains a good approxima- 
tion. 



form of the A cosmology (see [53, [47|). as expected. 
Obviously, for the Einstein-de Sitter model (Q m = 1, 
7 = 0) we have fl\ = and all coefficients vanish 
n c = n t = no = 0, so that eq. (|45[) boils down to 
A = g x (0) = 1/2, seeeq. (|43l) . 

In this framework, the density contrast at the virial- 
ization epoch is given by: 



Pr, 



Pr, 



A 3 \a t 



(49) 



where p„ 



ture and p m ^ c 



is the matter density in the virialized struc- 
is the background matter densit y a t the 
same epoch. Following the notations of J5J, IZ3, l74j . 
we again provide in appendix C a fitting formula for 
A V i r (within a physical range of cosmological parame- 
ters: < 7 < 0.01). Notice that the Einstein-de Sitter 
value for A„j r is precisely 18-7T 2 , and was factorized in 
CI 4 . 
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FIG. 3: Upper Panel: The deviation (1 - A/Aa)% of the col- 
lapse factor for various vacuum models with respect to the A 
solution. Bottom Panel: The density contrast at the virializa- 
tion, A v i r , as a function of redshift. The lines represent the 
following cosmological models: (a) Argh (solid), (b) Argci 
(dotted line) and (c) Argc2 (dashed line). The points repre- 
sent the traditional A cosmology. 



4 This value ensues from the Einstein de-Sitter value of £ = 
(37r/4) 2 multiplied by 2 s . Indeed, the sphere contracts a factor of 
2 from the turnaround point to virialization, and the background 



TABLE I: Numerical results. The 1 st column indicates the 
vacuum model used. Between column two and four, we 
present the main properties of the spherical collapse model 
[A„j r (z c ), (,{zt)], assuming that galaxy clusters have collapsed 
prior to the present time z c ~ (zt — 0.75). In columns four 
and five, we give the same quantities but considering that 
clusters have formed (collapsed) prior to epoch of z c ~ 1.6 
(zt — 3.2), in which the most distant cluster has been found. 



Model 


A„; r (0) 


C(0.75) 


A„; r (1.6) 


C(3.2) 


A 


348.6 


6.80 


190 


5.64 


Argh 


339 


6.79 


184.3 


5.62 


Argci 


317.6 


6.66 


168.6 


5.40 


ArgC2 


368.8 


6.62 


195.2 


5.38 



4.2. Clustered Vacuum Energy 

In this section we consider a scenario in which the vac- 
uum energy density on the scale of galaxy clusters is clus- 
tered: pAs{R) = A s (R)/8ttG. In such a scenario it could 
be possible, on non-linear scales, to have an interaction 
between dark matter and dark energy with a different 
7 than the background value. In the overdensity rest 
frame and in the homogeneous case (described before), 
the dark energy component flows progressively out of the 
overdensity [58., 62], and hence energy conservation can- 
not be applied (especially for large values of 7's) in or- 
der to determine the collapse factor A (along with the 
virial theorem) . To simplify the inhomogeneous case for- 
malism, we consider the extreme situation in which the 
vacuum energy fully clusters along with the dark matter, 
avoiding energy non-conservation which was examined in 

Within this framework, we also assume that the gen- 
eral functional form that describes the behavior of the 
vacuum energy density inside the spherical perturbation 
obeys a similar equation as that of eq. (f2TJ) : 



A 3 (R) =A q + 8ttG 



or else, 
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(50) 



(51) 



where the slope 7 S is not necessarily equal to the back- 
ground 7. It is worth noting that if we leave 7 S to take 
also negative values then there is a critical radius i?*, 
of the spherical overdensity, in which, for R < R+, the 
inhomogeneous vacuum energy density becomes negative 
A s < 0. Thus, the extra positive pressure (P\ s > 0) in- 
side the spherical perturbation implies compression (sim- 
ilar to gravity) rather than tension (while the opposite 
is true for R > i?*). With the aid eq. (|5T|) . this critical 
radius is given by 



scale factor increases a c /at = 2 2 ' 3 , thus (2 X 2 2 ' 3 Y 



:2 5 . 



i?,+ — R t 



1- 



"o(l -7s) 
27s 



'3TT^ 



(52) 



In the previous equation, we borrowed the definition of 

n given in (J48]), where p m& ,t = p m ,t( = p m oa t { 1 ~ 7) C is 
used, in which 7 (not 7 S ) is involved because it refers to 
the evolution of the background matter density. 

On the other hand, if 7 S > then A s > (or Pa s < 0) 
for all values of < y < 1. In this paper we use 2 
different versions of A s , namely 7 S = 7 = 0.002 (hereafter 
Argci) an d 7s = —0.002 (hereafter A^gc2)- Of course 
for the basic cosmological functions [il m (a), E(a) and 
D(a)] which enter in this vacuum pattern, we utilize the 
background 7 = 0.002. 

Inserting eq. tfSTj) into eq. (j37|) , we obtain 
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H t fl r 



(1 ~ 3 7 .)C 

2-37 s 



1 



is)y 



1st 

1-7, 



(53) 

In contrast with the homogeneous case, the novelty of the 
current approach is that the above differential equation 
can be solved analytically. Indeed, due to the fact that 
the differential eq. (|53[) is a function only of y we can 
perform easily the integration 



y 2 = H 2 t n m , t [P{y,Q + C] 
where C is the integration constant and 



P(y,0 = 



C 



7,C 



II 



(54) 



(55) 



(l- 7s ) 2 / 1 - 3 ^ \ i- 7s , 

Using now eq. (|40p we can provide the basic differential 
equation for the evolution of the overdensity perturba- 
tions 

(56) 



d.r 



P{yX) + c 



a;-i+37 4. rx 2 I(x) 



where the boundary conditions, (dy/dx) — and y — 
1 at x = 1, imply that C = — P(l,£). Therefore, the 
general integral equation which governs the behavior of 
the density contrast £ at the turn-around epoch, for the 
RG vacuum models is 



dy 



dx 



^P{yX) + C Jo ^x-i+37 + rx 2 I(x) 



(57) 



Of course, the usual ACDM cosmology is fully recovered 
from this model in the limit 7 = j s =0. 

Similarly, as in section 4.1, we again provide a useful 
fitting formula for £ as well as for A„j r , as a function 
of the cosmological parameters (for more details see ap- 
pendix C). 

Now with the aid of eq. (j5T))) . we can integrate eg . ([2!))) 
in order to derive the potential energy associated with 
the vacuum energy inside the spherical overdensity. In 
particular, we show here that U\ can be written as a 
sum of three components that contribute to the local 
dynamics 



Ua = _M^ R i + 4nG l?M£^ R i 

5(1-7 S ) 
J s Mp mStt 



10 
- 4ttG 



(l- 7s )(2 + 3 7s ) J Rr 3(1_7) 



ir 1+37 % (58) 



where p ms ,t — (p m ,t (see section 4). In this case, the 
algebraic equation which defines the collapse factor is 
found from the combination of equations (|42|) . (f28|) and 
(1551) as follows: 



g2(7,)[«o-/(7 s )]A 3 -A(n ,7,)A+.g(7 s )A 373 +2 gi (7 s ) - 0, 

(59) 
where uq is defined in equation (|48p, with 

,/ n 2 7s / s 107 s <72(7s) 

/(7s) = -, , 9V7s) = 



1-7* 



and 



A(n ,7 s )=2 + n -/(7 s ) + 



(l- 7s )(2-3 7s ) 



12 (is) 



Finally, solving eas. ([ST|) [or using ea lCl] and ([59^1 . we 
can estimate the density contrast at virialization from 
eq.gS). 



5. COMPARISON AMONG DIFFERENT TYPES 
OF VACUUM 

In this section, we investigate in more detail and com- 
pare the spherical collapse model, using different versions 
of the A(t) model (homogeneous or clustered), with that 
of the traditional A cosmology (see appendix A). In fi- 
gure 2 we present the density contrast at the turnaround 
point, Cj as a function of the turnaround redshift, Zt, for 
the constant vacuum A (solid points) , homogeneous vac- 
uum Arch (solid line), clustered vacuum Argci (dot- 
ted line) and clustered vacuum Arqc2 (dashed line). It 
is obvious that the A and A^qh models are almost in- 
distinguishable, while C, appears to be somewhat lower in 
the inhomogeneous (clustered) case. Indeed, we find that 
Cci/Ch ~ 7r/3 at large redshifts. 

Solving now eq. (|4"5|) and eq. (f5l?]) we calculate the col- 
lapse factor and we find that it lies, in general, in the 
range 0.46 < A < 0.52, in agreement with previous stud- 
ies ( H, m,"jH|, jg3|, (El, H, [11). In the upper 
panel of figure 2 we plot the deviation, (1 — A/Aa)%, of 
the collapse factors, A(z c ), for the current vacuum mod- 
els with respect to the A solution, Aa(z c ). It becomes 
evident that the size of the cosmic structures which are 
produced in the Arch model (solid line) is remarkably 
close to that predicted by the usual A cosmology, and 
therefore the impact of the vacuum energy on the spher- 
ical collapse is very small in the homogeneous case. This 
was to be expected. 

On the other hand, when considering the effect of the 
clustered vacuum energy, the largest positive deviation of 
the collapse factor occurs for the ArgC2 model (dotted 
line), which implies that this model produces more bound 
systems than the concordance ACDM model. There- 
fore, within this vacuum pattern the corresponding cos- 
mic structures should be located in larger density envi- 
ronments. The opposite situation holds for the Arqci 
(dashed line) model due to its negative deviations. 
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In figure 3 we plot the evolution of the density contrast 
at virialization. At very large redshifts, it tends to the 
Einstein-de Sitter value (A^ ~ 187r 2 ), as it should. 

In Table 1 we list, for the case of a cluster forming at 
z c ~ or at z c ~ 1.6, the following: (a) the cosmological 
models and the value of the turn around redshift; and 
(b) the virial density A V i r at the collapse time, as well as 
the density excess of the matter density in the spherical 
overdensity, £, at the turn around time. 

We also verify that the density contrast decreases with 
the formation (virialization) redshift z c . The A V i r differ- 
ences among the different vacuum models enter through 
£l m (z) [see eq. (|C3[) ] as well as through the assumption 
about the behavior of the vacuum inside the spherical 
overdensity (homogeneous or not). This feature points 
that perhaps the density contrast at virialization can be 
used as an effective cosmological tool. 

Although, we will investigate in detail such a possibil- 
ity in a forthcoming paper, we present here an idea of 
how to use observational data to estimate A vir . Using 
existing catalogs of clusters of galaxies, one should select 
those clusters which appear to have a quasi-spherical pro- 
jected shape 5 , which is the expected shape of a virialized 
cosmic structure, and then derive their virial radius and 
One can then easily calculate the observational 



mass. 



value of A. v i r of a cluster at a redshift z, from: 



A., 



3M„ 



47rp crit ,o^m,o(l + z) 3 r. 



3 

vir 



(60) 



and compare it with the model expectations. Now, the 
cluster virial radius can be calculated from the projected 
separations of the N m galaxy members according to (eg., 



N m {N m - 1) 



2 £?=i E;=! +1 [^tan(<^)]- 



(61) 



where di, is the luminosity distance of the group and 66ij 
is the angular (i,j)-pair separation. Using the observed 
cluster velocity dispersion, a v , and rw one can estimate 
the cluster's virial mass using the virial theorem, accord- 
ing to: 



M„ 



3a, 



G 



5G 



(62) 



The second A-based term is negligible, ~ 4.7 x 
10 11 fl\r^ ir Mq, and therefore it does not affect signif- 
icantly the mass estimates of clusters of galaxies. Of 
course, this approach is of a statistical nature, since there 
are various observational systematics that enter in the in- 
dividual cluster determination of A V i r , as well as cosmic 



5 Not all apparently spherical clusters are truly spherical since 
elongated clusters with their major axis oriented at small an- 
gles along the line of sight will appear spherical in projection, a 
fact which is a further source of noise. 



variance, which however can be minimized if one averages 
over a suitable and relatively large sample of clusters at 
each redshift interval. 
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FIG. 4: The present time virial density for the clustered vac- 
uum energy model as a function of 7 S (open points). The 
inner and outer dashed lines correspond to the 1 and 2a, 
observationally determined, virial density limits respectively 
(202 < A v i r < 494), based on a subsample of the 2MASS 
High Density Contrast group catalog. 

We have applied this method olog y to the 2MASS High 
Density Contrast group catalog [7f|, which is a low- 2 cat- 
alog based on the 2 micron infrared whole-sky survey and 
which was constructed by a friends-of-friends algorithm 
(eg. [771] ) such that the groups correspond to an overden- 
sity Sp/p > 80. This catalog was carefully constructed, 
with respect to other catalogs, and it is less prone to 
projection, interloper contamination and contamination 
by the large-scale structures from which galaxies are ac- 
creted to the groups (see [78[ for a relevant discussion). 
We selected only those groups with projected axial ra- 
tio > 0.8 and with at least 16 galaxy members (in or- 
der to have a relatively accurate determination of their 
shape, velocity dispersion and thus M V i r ) and we are 
thus left with 7 clusters at (z) ~ 0.015. We clip the 
lower and higher A V i r outliers, since we do expect sys- 
tematic effects to be present, and we derive a mean value 
of (A v i r ) = 348 and a standard deviation of the distribu- 
tion of ±73 (if we use N m > 20 we are left with 6 clusters 
with (A V i r ) = 329±69). Although, we have derived these 
A V i r values using the concordance ACDM cosmological 
model to estimate cLl , there would be no appreciable dif- 
ference had we used any of the other models, presented 
in section 2.1 (because of the very small value of 7 and 
of the very low redshift of the sample). Inspecting Table 
1 it is evident that the previously derived observational 
values are in good agreement with the theoretical expec- 
tations although with the present level of uncertainty we 
cannot distinguish among the models. However, in the 
case of the clustered vacuum energy model we can put 
some limits on the value of 7 S even with present level of 
accuracy. As an example, we compare in Figure 4, the 
predicted virial density, A vir , of the model at the present 
time (open points), with the observationally derived 2<r 
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range of values, based on the previously discussed sub- 
sample of the 2MASS High Density Contrast group cat- 
alog, and find a consistency for —0.009^ 7 S <; 0.012. In 
the future we plan to further investigate the model pre- 
dictions, using a larger number of clusters spanning also 
a range of different redshifts, in an attempt to put strin- 
gent constraints on the value of 7 S . 



and 
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il\ + YL m a 



(A2) 



(A3) 



6. CONCLUSIONS 

In this paper we have studied analytically and numer- 
ically the spherical collapse model in the case of a time 
varying vacuum, with A(H) = no + ri2 H 2 , for a spatially 
flat FLRW geometry. We find that the amplitude and the 
shape of the virial density contrast is affected by the con- 
sidered status of the vacuum energy model (homogeneous 
or clustered) . We verify that in the case where the distri- 
bution of the vacuum energy is clustered the structures 
produced are more concentrated (under specific condi- 
tions) with respect to the homogeneous dark energy case. 
Finally, by comparing the predicted virial density con- 
trast at the present epoch with a preliminary analysis of a 
suitable subsample of the 2MASS High Density Contrast 
group catalog (at a mean redshift of (z) ~ 0.015), we 
find that the inhomogeneous vacuum energy models can 
be accommodated, at a 2a level, if the vacuum clustering 
parameter is within the range: —0.009^ 7 S <^ 0.012. The 
latter result points to the direction that perhaps the A V i r 
parameter, once estimated accurately from observations, 
could be used in order to determine the internal physical 
properties of the vacuum energy. 
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• The spherical model: The basic set of equations 
here is: 



# t 2 ft ro , t [aT 1 + rx 2 I{x)] 



y 



H t vt m ,t f C 



— - 2ryl(x] 



2 \y 2 
where I(x) = 1 and 

"A,i 3 

r = -— = r a t 



Q 



(A4) 



(A5) 



(A6) 



rn , / 



where r$ = VL^/Q, m . Therefore, the general integral 
equation which governs the behavior of the density 
contrast £ at the turn around epoch is 



dy 



lnU/TT7 + V?) 2 / 3 



o y/P(y,C)-P(l,0 ^ 

where 

P(yX) = - + ry 2 . 
y 



(A7) 



(A8) 



Note, that the time needed for a spherical shell to 
collapse is twice the turn-around time, tf ~ 2t t . 
This implies that: 



sinh x ( \Jr af. ) ~ 2sinh x 



r at 



(A9) 



Appendix A: The Concordance A cosmology 

In this appendix we would like to give the reader the 
opportunity to appreciate the fact that the A^g model 
can be viewed as an extension of the concordance A cos- 
mology. In particular, the basic cosmological equations 
in the A/jg model reduce to those of the A cosmology for 
7 = 0. Below, we present the main quantities of the A 
cosmology: 

• Global Dynamics: The basic cosmological equa- 
tions (see section 2.1) take the following forms: 



(Al) 



H(t) = V^a H o coth 



• Virial Theorem: The virial theorem becomes: 

2T + U G - 2U A = . (A10) 

Using now also the energy conservation at the turn 
around and at the virial time we derive the follow- 
ing relations: 



1 



U G ,c + 2t/ A , c = Ua.t + U A . t 



(All) 



2n A 3 - (2 + n )A + l = (A12) 

where A = R c /Rt is the collapse factor and 



no 



2f7 A a 3 



(A13) 
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Appendix B: Roots of a cubic polynomial 



where 



We remind the reader of some basic elements of Alge- 
bra pertinent to our analysis. Given a cubic equation: 
A 3 + a x A 2 + a 2 A + a 3 =0, let D be the discriminant: 



V 



2 2 

a 1 a 2 



4a| 



and 



4a?a 3 



27 



27a\ + 18aia 2 a 3 (Bl) 



3V3X> 



X\ = -a x + -a x a 2 — a 3 



X2 



If 2? > 0, all roots are real (irreducible case). In that case 
Ai, A 2 and A 3 can be written: 



Ol 

' 3 



2r l/3 



-cos 



(p - 1)tt 



H = 1,2,3 



^ r (B2) 

where r = y^i + ^2 an< ^ ^ = cos 1 (a;i/r). 

Now, we are ready to derive analytically the exact 
roots of the basic cubic equation (|45j) having polyno- 
mial parameters: a x = 0, a 2 = — (2 + n t )/q 2 n c and 
a 3 = 2qi/q2n c . Then the discriminant becomes: 



V(n t ,n c ) 



(2 + n t f - 21q\q 2 n c 



q\n% 



(B3) 



Of course, in order to obtain physically acceptance solu- 
tions we need to take n t ,n c > 0, which gives V{n t ,n c ) > 
0. Therefore, all roots of the cubic equation are real (ir- 
reducible case) but one of them < A 3 < 1 corresponds 
to expanding shells. It is obvious that for n t ,n c — > 0, 
the above solution tends to the Einstein-de Sitter case 
(A3 — > 0.50), as it should. 



Appendix C: Fitting formulae 

We provide here accurate fitting formulae for the den- 
sity contrast at the turn around redshift and at the epoch 
of virialization, and which do not have a simple fully an- 
alytical form. These are: 



3tt 



Vi + A. 



n 



-a>l+W2fi m ,t— U)3W(<h) 



(CI) 



w(a) = — 1 — 



7«' 



37 



a 3 f + r 



(C2) 



and 



A vir (a)~ 187r 2 [l + e6 6 (a)] 
where 9(a) = Cl^(a) - 1. 



(C3) 



• Homogeneous Vacuum: In this case we have A s 
0, (wi,w 2 ,W3) = (0.79,0.26,0.06) and 

e = 0.40 - 25 7 + 500 7 2 , 6 = 0.94 + 507. 

• Clustered Vacuum: Here we find 



A = f -24.25 7s + 2125 7s 2 < 7s < 0.01 ,_ 

8 I 29.75 7s + 2375 7s 2 - 0.01 < 7, < 



(W1,W2,W3) 



(0.86,0.36,0) 0<7 S <0.01 
(0.74,0.16,0) -0.01<7 S <0 



(C5) 



I 0.94 + 145 7s + 4.75 x 10 4 7 2 - 0.002 < 7s < 0.01 

[ 0.94 + 557, - 0.01 < 7, < -0.002 

(C6) 
and 

f 0.40 - 65 7;s - 1.25 x 10 4 7 2 - 0.002 < 7s < 0.01 
[ 0.31-86.257 s -0.01 <7 S < -0.002. 

(C7) 
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